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We investigate the spin structure of the pion within the framework of the nonlocal chiral quark 
model from the instanton vacuum. We first evaluate the tensor form factors of the pion for the 
first and second moment (n = 1,2) and compare it with the lattice data. Combining the tensor 
form factor of the pion with the electromagnetic one, we determine the impact-parameter dependent 
probability density of transversely polarized quarks inside the pion. It turns out that the present 
numerical results for the tensor form factor as well as those for the probability density are in good 
agreement with the lattice data. We also discuss the distortion of the spatial distribution of the 
quarks in the transverse plane inside the pion. 
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I. INTRODUCTION 

The transversity of hadrons has been one of the most important issues well over decades (see a recent review [I]), 
since it allows one to get access to their spin structures. It is pertinent to the tensor current of hadrons and is very 
difficult to be measured experimentally, because there is no direct probe to measure it. Only very recently, it was 
suggested that the transverse spin asymmetry Att in Drell-Yan processes in pp reactions [2H5] as well as the azimuthal 
single spin asymmetry in semi-inclusive deep inelastic scattering (SIDIS) [B] can be used to obtain information on the 
transversity of the nucleon. Though it is even more difficult to measure the transversity of the pion experimentally, 
it is still of great significance to understand it, since it provides the internal spin structure due to quarks, i.e. it 
accommodates a novel concept called a hadron tomography. The first result of the pion transversity on lattice has 
been reported by the QCDSF/UKQCD Collaborations 7 . They also presented the probability density of the polarized 
quarks inside the pion, combining the electromagnetic form factor of the pion [8] with its tensor form factor. It was 
shown in Ref. |7] that when the quarks are transversely polarized, their spatial distribution is strongly distorted. This 
first result in the lattice QCD has triggered several subsequent theoretical works In Ref. [TT], the tensor form 

factors of the pion have been studied within the local and nonlocal Nambu- Jona-Lasinio (N JL) model [TT] , a direct 
comparison with the lattice results being emphasized. In doing so, they employed a larger value of the pion mass, i.e. 
TO^ = 600 MeV such that the results can be confronted with the lattice data. They also considered the case of the 
chiral limit. 

In the present work, we first want to investigate the pion tensor form factor in the space-like momentum transfer 
region (0 < < 1 GeV), based on the low-energy effective chiral action (E^A) from the instanton vacuum jl2j . 
Combining the result of the tensor form factor with the electromagnetic one of the pion which was already studied in 
Ref. jl3j within the same framework, we then derive the probability density of transversely polarized quarks inside the 
pion. Since the instanton vacuum realizes the spontaneous chiral symmetry breaking (SxSB) naturally via quark zero 
modes, it may provide a good framework to study properties of the pion such as the electromagnetic and tensor form 
factors. Moreover, an important merit of this approach lies in the fact that there are only two parameters, that is, 
the average (anti)instanton size p ~ 1/3 fm and average inter-instanton distance _R 1 fni. The normalization point 
is given by the average size of instantons and is approximately equal to « 0.6 GeV. The values of the p and R 
were estimated many years ago phenomenologically in Ref. [14] as well as theoretically in Ref. [TSHTT] . The instanton 
framework has been proved to be reliable in reproducing experimental data especially for the meson sector, such as the 
meson distribution amplitudes jl8H20| . semileptonic decays [21] . and etc. Furthermore, this approach was supported 
by various lattice simulations of the QCD vacuum ^22n24] . The quark propagator from the instanton vacuum | 15j 
is in a remarkable agreement with lattice calculations [551 [2S]- Finally the nonlocal chiral quark model from the 
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instanton vacuum has a practical virtue, since it does not have any adjustable parameter once the above-mentioned 
two parameters p and R are determined. 

We organize the present work as follows: In Section II, we briefly explain the definitions of the probability densities 
of the transversely polarized quarks, which arc expressed in terms of generalized form factors of the pion. In Section 
III, we show how to calculate the tensor form factors within the nonlocal chiral quark model from the instanton 
vacuum. In Section IV, the numerical results are discussed and compared with those in lattice QCD. The final 
Section is devoted to summarize the present work, to draw conclusions, and to give outlook. 



II. GENERALIZED FORM FACTORS OF THE PION 



In this Section, we define the probability density of transversely polarized quarks inside the pion. For definiteness, 
we choose the positively charged pion tt"*" from now on. The probability density of transversely polarized quarks for 
the n th moment of the probability density is given as 



dx : 



(1) 



where denote the impact parameter that measures the distance from the center of momentum of the pion to the 
quark in the transverse plane to its motion. The sj_ stands for the fixed transverse spin of the quark. For simplicity, 
we choose the z direction for the quark longitudinal momentum. The x indicates the momentum fraction possessed 
by the quark inside the pion. The An{){b\) and i3„o(^D are called the generalized form factors (GFFs). In fact, 
the GFFs are just the moments of the generalized parton distributions (GPDs) for the unpolarized and transversely 
polarized pions, respectively: 



dxx''-^H{xX = Q.bi)^ ArM{b\), 



j dxx''-'^E{x,C^O,bl) 



Bnoibl). 



(2) 



For the first moment, the GFFs Aiq and Biq are identified with the electromagnetic and tensor form factors of the 
pion, respectively. Previously, we have studied v4io((7^) the momentum space within the nonlocal chiral quark 
model (NL^QM) from the instanton vacuum IB], resulting in a good agreement with the experimental data. Hence, 
we can readily calculate AiQ{b±), using the results of Ref. PB]. Thus, we will concentrate on calculating the tensor 
form factors i?io.20 within the same framework, and they can be written in a general form as follows: 



q^'---q^'P^'+'---p^"-'B,,,iQ^) 



t—Gvcn 



(3) 



where pi and pf stand for the initial and final on-shell momenta of the pion, respectively. We also use notations 
p = {pf -\-pi)/2 and q ~ pf ^ Pi. The tensor operator also can be given as: 



O 



AS q''a^"'{iD^'')---{iD''--')q 



(4) 



The A and S denote the anti-symmetrization in (/i, J^) and symmetrization in (i^, ■ ■ • , fin-i) with the trace terms 
subtracted in all the indices. Taking into account Eqs. ([s]) and Q, we can define the tensor form factors Biq and B20 
of the pion in momentum space as the matrix elements of the tensor current, using the auxiliary- vector method as in 
Ref. Ell: 



(7r+(p/)|gt(0)a,bg(0)|7r+(K)> = [ip^ ■ a)ipf ■ b) - (p, • b){pf ■ a)] 



{n+{pf)\q\0)aabiS ■ a)g(0)|^+(K)) = {{p ■ a)[(p, • a)ipf ■ b) - [p, ■ b){pf ■ a)]} 



B2o(Q') 



(5) 



where the vectors satisfy the conditions, i.e. — 0, a ■ b — and b^ ^ 0, and we have used a notation Oah = cr^,ya^6''. 
Due to this auxiliary-vector method, one can eliminate the trace-term subtractions. We also introduce a notation 
= (iD ^ — iD^)/2^ where indicates the SU(A^c) covariant derivative. Since we are interested in the spatial 
distribution of the transversely polarized quark inside the pion, we need to consider the Fourier transform of the form 
factors: 



1 /■ 1 

F{bl) = j d^^e-^'^-^^Fiql) = — QdQMbQ)F{Q^), 



(6) 
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where F = {Aio,Bio) designates the generic pion form factor. The magnitudes of the transverse momentum and 
impact parameter are expressed as \q±\ = Q and |b_L| = b. Similarly, the Fourier transform of the derivative of the 
GFF with respect to 6^ can be evaluated as: 



dF{bl) _ , 1 



Q^dQJi{bQ)F{Q^). 



(7) 



The Jo and Ji in Eqs.(6[7l denote the Bessel functions of order and 1, respectively. According to the definitions for 
the relevant vectors and the probability density in Eq. ([T]) reads as follows: 



Pi {b±,Sr, = ±1) = - 



Aio(6') T 



bsinO 



(8) 



where the spin of the quark inside the pion is quantized along the x axis, s± — (±1,0). 



III. NONLOCAL CHIRAL QUARK MODEL FROM THE INSTANTON VACUUM 

We now briefly explain the NLxQM from the instanton vacuum ,28J and derive the GFFs of the pion. Considering 
first the dilute instanton liquid, characterized by the average (anti)instanton size p ~ 1/3 fm and average inter- 
instanton distance R ^ 1 fm with the small packing parameter np'^/R'^ « 0.1, we are able to average the fermionic 
determinant over collective coordinates of instantons with fermionic quasi-particles, i.e. the constituent quarks intro- 
duced. The averaged determinant is reduced to the light-quark partition function that can be given as a functional of 
the tensor field in the present case. Having bosonized and integrated it over the quark fields, we obtain the following 
effective nonlocal chiral action in the large Nc limit in Euclidean space: 



5eff [m, tt] = -Sp In Up + im + i^/M{id)W- (</>) ^/M{id) ■ 



■T 



(9) 



where m, tt, and Sp indicate the current quark mass, the Nambu-Goldstone (NG) boson field, and the functional trace 
over all relevant spaces, respectively. In the numerical calculations, we will choose m ~ to„ ~ rud ~ 5 MeV, taking 
into account isospin symmetry. The M{id) stands for the momentum-dependent effective quark mass, generated from 
the nontrivial quark-(anti)instanton interactions [12J. Although its analytical form is in general given by the modified 
Bessel functions, we will make use of its parametrization for numerical convenience: 



M{id) = A/o 



2 + fP-d^ 



(10) 



where Mq indicates the constituent quark mass, which can be determined self-consistently by solving the gap equation 
of the present framework, resulting in Mq = 350 MeV [T^]. The NG boson field is represented in a nonlinear form 
as ESI: 



[/'^^(</)) = exp 



»75(t ■ </>) 



2Fi 



(11) 



where (j)"" is the SU(2) multiplet, defined as 



T ■ (p = 



V2 



(12) 



The F^ denotes the weak-decay constant for NG bosons, whose empirical value is 93.2 MeV for the pion for instance. 
The last term in Eq. ([9| denotes cr • T = cr^iyT)ii/, where ct^^ = ^[7^, 7i/]/2 and represents the external tensor 
source field. 

The three-point correlation function in Eq. ([5| can be easily calculated by a functional differentiation with respect 
to the pion and external tensor fields, which leads to the following two terms for the Bio{Q^): 



ST Stt'^Stt'' 



T=0 



p2 



Sp 



Af 75 r " V A/ — VMj^tWM 



2Fi 



Sp 



^ fj.iy 



, (13) 
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where we have introduced the shorthand notations M = M{id) and ip = ip + im + iM(id) = ip + iM{id). The 
trace over the isospin space yields tr^ [r^r''] = 25°'^ . One can also do for the i?2o(Q^) similarly. Having performed the 
functional trace and the trace over the color space, we arrive at the matrix elements for the -B(io,20) {Q^)i corresponding 
to Eq. ([5]), as follows: 



(^+(p/)|gt(0)a,b<;(0)|7r+fe)) = - 



2N^ f d^k 



r+(p/)|gt(0)a,,(iS.a)q(0)|7r+fe)) 



2Nr. 



J (27r) 



4^1-7 



(A) 



Fl J (2^) 



-Tr 



(B) 

tPa iPb 2 



(A) 

iPb ^Pc 2 



(14) 



(B) 



The corresponding Feynman diagrams to the two terms (A) and (B) in the right-hand side of Eq. ( 14 1 are depicted 
in Fig. [1] respectively. The relevant momenta are also defined as: 





FIG. 1: Feynman diagrams for Bno in the NLxQM. We assign the initial and final pion momenta as pi and p/, respectively, 
while the momentum transfer as q. We also define the loop momenta aska = k— — %,kh = k+ ^ — %, and fcc = fc + ^ + f . 



ka — ^ ~ ~2 ~ 2' ^'^ ~ ^ ^ ~2 ^ 2' ^'^ ~ ^ ^ ~2 ^ 2' 



(15) 



In order to evaluate the matrix element, we define the initial and final pion momenta in the Breit (brick- wall) frame 
in Euclidean space as done in Ref. |13) : 



-|,0,0,i^^+m2 j , ( |,0,0,z^^+m2 j , g=(Q,0,0,0), 



(16) 



We also have chosen the auxiliary vectors for definiteness as a = (0, l,0,i) and b = (1, 0, 1^ 0), which satisfy the 
conditions mentioned in Section II. The denominators become pa,b.c = j^ a,b,c + iMahc in Eq. (14). The momentum- 
dependent effective quark mass Ma^b,c can be also defined by using Eqs. (10) and (15). 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



We first discuss the numerical results of the tensor form factors of the pion. Figure [2] draws the numerical results 
for the electromagnetic form factor of the pion Aiq in the left panel and its tensor form factor Biq in the right panel 
as functions of in the range 

of < Q2 < 1 GeV^ However, we want to mention that there is a caveat. Since we 
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FIG. 2: The numerical results for the electromagnetic form factor Aio{Q ) are presented in the left panel and those for the 
tensor form factors -Bio(Q^) and B2o{Q^) in the right panel. The solid curves depict the numerical results, whereas their 
parametrizations given respectively in Eqs. (171 and (211 are denoted by the squares. We also show the parametrized lattice 
data [7] for the tensor form factors, designated by the triangles in the right panel. 



need the results of the form factors in principle up to infinite in order to perform the Fourier transform given in 
Eq. (|6|, we will use the parametrized one, as we will discuss later in the context of the lattice data. 

The numerical results for Aio{Q^) are taken from Ref. [T3j. Though we have already discussed those for the 
electromagnetic form factor in detail in Ref. |13j . we want to recapitulate them in the context of the lattice data. It 
is well known that the electromagnetic form factor can be parametrized by a monopole form 

F.(Q^)^A,o(Q^)^^^^. (17) 

The monopole mass M was known to be M = (0.714 ± 0.004) GeV, based on the experimental data [53H3T]. On the 
other hand, the lattice calculation yields M = (0.727 ± 0.016) GeV with a linear chiral extrapolation to the physical 
pion mass taken into account 8J. The present result leads to M = 0.738 GeV, which indicates that that of the pion 
electromagnetic form factor is in good agreement with the lattice data. We also obtain the squared charge radius of 
the pion (r^) = 0.456 fm^, while in the lattice QCD it was evaluated to be (r^) = (0.441 ± 0.019) fm^. Conside ring the 
uncertainty of the lattice data, the present result is in remarkable agreement with them. In the left panel of Fig. [2j we 
show the numerical result (solid curve) [13' and its monopole parametrization (sqaure) of Aiq{Q^), using the values 
mentioned above and Eq. ( |17t . 

In the right panel of Fig^M we draw the numerical results for the tensor form factors -Bio(Q^) and i?2o(Q^) (solid 
curve). In order to compare the present results with the lattice data, it is crucial to consider the evolution of the 
scale [T1[TT1[32], since the tensor current is not the conserved one. The tensor form factor is evolved at the leading 
order (LO) by the following equation [TJ [TT] 



-B„o(Q',m) = S„o(Q',Mo) 



a{po) 



(18) 



where we have used the anomalous dimensions 71 — 8/3 and 72 — 8, and /3o — llA'c/3 — 2Nf/3 {Nc — 3 and Nf = 2 
in the present case). Thus, the powers in the LO evolution equation are given as 4/29 and 12/29 respectively for 
n = 1 and n = 2, which indicate that the dependence of the tensor charge on the normalization point turns out to 
be rather weak. Note that the anomalous dimension is simply the same as that for the nucleon tensor charge |33) . 
We also take Aqcd = 0.248 GeV which was also used in evolving the nucleon tensor charges and tensor anomalous 
magnetic moments |331I3S]. Since the normalization point of the present model is around 0.6 GeV, while the lattice 
calculation was carried out a,t fi — 2 GeV, the scale factors turn out to be 

-Bio(Q',M = 2GeV) = 0.89Blo(Q^/io = 0.6 GeV), B2o(Q',m = 2 GeV) = 0.70 Ban (Q',Aio = 0.6 GeV). (19) 

Gonsidering these scalings, we obtain i3io(0) = 0.216 and i?2o(0) = 0.032. In the lattice calculation [7 , the tensor 
charge of the pion for n = 1 with the linear chiral extrapolation to the physical pion mass in was estimated to 
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FIG. 3: -Bio(0)/m^ (left) and -B2o(0)/m^ (right) as functions of m^. The numerical results for = and 140 MeV are 
given by the square and circle. The lattice data (triangle) are taken from Ref. [2 . The shaded areas for m,r = (0 ~ 140) MeV 
represent the lattice fits. 



be about i?io(0) = 0.216, which is almost identical to the present result. As for n = 2, the lattice data estimated 
about B2o{0) = 0.039, which is about 20% larger than the present one, but is still comparable. We want to mention 
that one could use larger current-quark masses in order to compare directly with the lattice data as done in Ref. jJT] . 
However, it is rather unreliable in the present framework: Firstly it is nontrivial to include the larger current quark 
mass |36H38j . Secondly, the present scheme is conceptually only valid when the current quark mass is small, at least 
up to the strange current quark mass. Thus, the NL^QM from the instanton vacuum is a rather restricted one, so 
that we will compare the present results with those of the lattice QCD with chiral extrapolation. 

In Fig. [3] we present the m.^- dependence of the pion tensor form factor scaled by the pion mass as a function of 
for ttItt = (square) and 140 MeV (circle). As shown in Fig. [sj the result in the chiral limit is slightly smaller than 
that with m^r = 140 MeV. As for the case with = 0, we take the current quark mass m — 0. The shaded bands 
represent the fits from the lattice calculation [T . 

A simple p-pole parametrization of GFFs was used in Ref. [7] to get the tensor form factor of the pion: 



BnoiQ^) = 5no(0) 



(20) 



In this parametrization, the lattice QCD simulation estimated the pole mass rrip^ — (0.756 ± 0.095) GeV and mp^ — 
(1.130±0.265) at = 140 MeV with chiral extrapolation. Considering the condition p > 1.5 for the regular behavior 
of the probability density at 6j_ — >■ |39] and following Ref. [7], we take pi = p2 = 1-6 as a trial. If this is the case. 



Eq. (20) gives us ra^^ 0.761 GeV and rup^ ~ 0.864 GeV to reproduce the present results, which is compatible with 

Taking into account these results, we can write the p-pole parametrized tensor form 



that of the lattice simulation, 
factor as follows: 



-Bio(Q') - 0.216 



1.6 X 0.7612 GeV^ 



-1.6 



S2o(Q') = 0.032 



1.6 X 0.8642 Gcv2 



-1.6 



(21) 



The result of this parametrized one in Eq. (21) is also depicted in the right panel of Fig. [2] (square) , and reproduces 
well the present numerical one. We also compare our result with the lattice one in the right panel of Fig. [2] As in the 
case of the electromagnetic form factors, the present results are in excellent agreement with the lattice data (triangle). 
Note that in the present framework we do not have any adjustable free parameter. 

We are now in a position to discuss the results of the probability densities of the transversely polarized quarks 
inside the pion, defined in Eq. ([T]) . In the upper-left panel of Fig. |4j we show the unpolarized probability density with 
the tensor form factor turned off. As expected, the quarks are distributed symmetrically on the bx-hy plane. On the 
other hand, if we switch on the tensor form factor, the spatial distribution of a transversely polarized quark inside 
the pion (tt"*") gets distorted as shown in the upper-right panel of Fig. [4] Its maximum value is also shifted to the 
by direction in comparison to that for the unpolarized one. Thus, it is interesting to examine the average transverse 
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-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



bx [fm] bx [fm] 




FIG. 4: In the upper panels, we show unpolarized (left) and polarized (right) probability densities, pi in Eq. (|T| as a function 
of and by for s±_ = (+1,0). Similarly, we show their three dimensional profiles in the lower panels. 



shift which is defined as [7l: 



where we have chosen s±_ = (+1, 0). Using Eq. (21 1 and m^r = 140 MeV, we obtain {hy) = 0.152 fm, which is almost 
the same as that of the lattice calculation {hy) — (0.151 ± 0.024) fm. This finite value of {hy) measures how much 
the polarized probability density is distorted in the transverse plane. If we take the spin quantized along the y axis, 
i.e. s_L = (0, +1), the result of the polarized probability density is similar but rotated by 90° clockwise. We also note 
that the present results are almost equivalent to those given by the lattice simulation ^7, . In the lower panel of Fig. |4j 
we show the three-dimensional profiles for the unpolarized (left) and transversely polarized (right) distributions as 
functions of h^ and by. One can obviously see that the maximum of the transversely polarized probability density is 
shifted and distorted. 

In Fig. [5] we draw the probability density as a function of by at bx ~ 0.2 fm, comparing it with that of the lattice 
calculation. As expected, the present result is almost identical to that of the lattice QCD. We summarize the main 
results of the present work in Table [l] 
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Q I I .... I .... I .... I .... I .... I .... I 

-0.3 -0.2 -0.1 0.1 0.2 0.3 

by [fm] 



FIG. 5: Comparison between the polarized probability densities from the present (square) and lattice (triangle) results for 
^ 0.2 fm. 

V. SUMMARY AND CONCLUSION 

In the present work, we have aimed at investigating the tensor form factors of the pion, Biq and B20, using the 
nonlocal chiral quark model from the instanton vacuum. Combining it with the electromagnetic form factor of the 
pion [13], we were able to evaluate the transversely polarized density of quarks inside the pion as functions of {bx, by) 
in comparison with the lattice simulation 7J. 

We first recapitulated the electromagnetic form factor of the pion computed previously in the context of the lattice 
calculation. We found that the monopole mass is M — 0.738 GeV which is in a good agreement with that from the 
lattice QCD M = 0.727 GeV as weh as the experimental data M = (0.714 ± 0.004) GeV. It indicates that the 
dependence of the electromagnetic form factor is well reproduced within the present work and compatible with the 
lattice results. We also presented the squared charge radius of the pion {r'^) — 0.456 fm^ which is again in very good 
agreement with the lattice result (r^) = (0.441 ± 0.019) fm^. 

We calculated the tensor form factor of the pion within the same framework as done in previous works. In 
order to compare the results with the lattice data, we evolved them from the fi = jj^^ — 600 MeV to the scale 
at which the lattice calculation was performed (/i = 2 GeV). We also carried out the p-pole parametrization as 
in the lattice QCD. The results for the tensor form factor were obtained as follows: The tensor form factors of 
the pion at = q (-B(io)(0), -B(2o)(0)) = (0.216, 0.032) and the pole mass mp^^^, = (0.762, 0.864) GeV. Being 
compared to the lattice results with chiral extrapolation, i.e. i3(io.20)(0) — (0.216 ±0.034, 0.039 ±0.099) and mp^j^ 2) = 
(0.756 ± 0.095, 1.130 ± 0.265) GeV, they were found to be almost identical and comparable to those of the lattice 
QCD. In particular, these results are remarkable, considering the fact that the present scheme does not contain any 
adjustable parameter. 

Having combined the results of the tensor form factor with those of the electromagnetic one, we obtained straight- 
forwardly the probability density of transversely polarized quarks inside the pion. It turned out that the spatial 
distribution of the quarks on the transverse plane were distorted, compared to that of the unpolarized quarks. More- 
over, the maximum value of the density is shifted to the by direction. In order to examine this shift, we also calculated 
the average value of by which turned out to be (by) = 0.152 fm. It is in an excellent agreement with the lattice result 
(by) = 0.151 fm. 

Finally, It is worth mentioning that it is also of great importance to study the spin structure of the kaon, since it 



rrin = 140 MeV 


Bio(O) nip, [GeV] (by) [fm] 


B2o(0) mp, [GeV] 


Present work 


0.216 0.762 0.152 


0.032 0.864 


Lattice QCD 21 


0.216 ± 0.034 0.756 ± 0.095 0.151 


0.039 ±0.099 1.130 ±0.265 



TABLE I: Numerical results for B„o{0), mp^, and (by) in comparison with those of the lattice simulation [7 with the chiral 
extrapolation to the physical pion mass — 140 MeV. 
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sheds light on the role of flavor SU(3) symmetry breaking inside the kaon. Related works are under progress and will 
appear elsewhere. 
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